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Abstract. The main objective of this article is to establish a new mecha- 
nism of the El Nino Southern Oscillation (ENSO), as a self-organizing and 
self-excitation system, with two highly coupled processes. The first is the os- 

Q dilation between the two metastable warm (El Nino phase) and cold events 
(La Nina phase), and the second is the spatiotemporal oscillation of the sea 
surface temperature (SST) field. The interplay between these two processes 
^v^j gives rises the climate variability associated with the ENSO, leads to both the 

random and deterministic features of the ENSO, and defines a new natural 

1 | feedback mechanism, which drives the sporadic oscillation of the ENSO. The 

new mechanism is rigorously derived using a dynamic transition theory devel- 
O | oped recently by the authors, which has also been successfully applied to a 

wide range of problems in nonlinear sciences. 

1. Introduction 

This article is part of a research program initiated recently by the authors on 
O ! dynamic transitions in geophysical fluid dynamics and climate dynamic. The main 

objective is to study the interannual low frequency variability of the atmospheric 
I and oceanic flows, associated with typical sources of the climate variability, includ- 

ing the wind-driven (horizontal) and the thermohaline (vertical) circulations (THC) 
^sO of the ocean, and the El Nino Southern Oscillation (ENSO). Their variability, in- 

dependently and interactively, may play a significant role in climate changes, past 
and future. 

ENSO is the known strongest interannual climate variability associated with 
strong atmosphere-ocean coupling, which has significant impacts on global climate. 
ENSO is in fact a phenomenon that warm events (El Nino phase) and clod events 
(La Nina phase) in the equatorial eastern Pacific SST anomaly, which are associated 
!* with persistent weakening or strengthening in the trade winds. 

. £^ It is convenient and understandable to employ simplified coupled dynamical 

models to investigate some essential behaviors of ENSO dynamics; see among many 

others [Ml E El [2H HU 13 HU EH1 [HI El 1231 [221 12B- 

An interesting current debate is whether ENSO is best modeled as a stochastic 
or chaotic system - linear and noise-forced, or nonlinear oscillatory and unstable 
system |19j ? It is obvious that a careful fundamental level examination of the 
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problem is crucial. The main objective of this article is to address this fundamental 
question. 

By establishing a rigorous mathematical theory on the formation of the Walker 
circulation over the tropics, we present in this article a new mechanism of the ENSO. 
Here we present a brief account of this new mechanism, and refer the readers to 
Section 5 for more precise description. 

We show that interannual variability of ENSO is the interplay of two oscilla- 
tion processes. The first is the oscillation between the metastable warm event 
(El Nino phase), normal event, and cold event (La Nina phase). We show that 
each metastable state has a basin of attraction, and the uncertainty of the ini- 
tial states between these basins of attraction gives rise of the oscillation between 
these metastable states. The second is the spatiotemporal oscillation of the sea 
surface temperature (SST) field, which is mainly caused by the solar heating and 
the oceanic upwelling near Peru. 

The interplay between these two processes give rises the inter-annual variability 
associated with the ENSO. On the one hand, the metastable states oscillation 
has direct influence on the SST, leading to 1) the intensification or weakening of 
the oceanic upwelling, and 2) the variation of the atmospheric Rayleigh number. 
Consequently, forcing the system to adjust to a different metastable phase, and to 
intensify/ weaken the corresponding event. On the other hand, the oscillation of 
the SST has a direct influence on the atmospheric Rayleigh number, leading to a 
different El Nio, La Nina and normal conditions. 

This interplay The above mechanism of ENSO as an interplay of the two pro- 
cesses defines a new natural feedback mechanism, which drives the sporadic oscilla- 
tion of the ENSO, and leads to both the random and deterministic features of the 
ENSO. The uncertainty is closely related to the fluctuations between the narrow 
basins of attractions of the corresponding metastable events, and the determinis- 
tic feature is represented by the deterministic modeling predicting the basins of 
attraction and the SST. Hence ENSO can be considered as a self-organizing and 
self-excitation system. 

The main technical method is the dynamical transition theory developed recently 
by the authors. The main philosophy of the dynamic transition theory is to search 
for the full set of transition states, giving a complete characterization on stability 
and transition. The set of transition states -physical "reality" - is represented by 
a local attractor. Following this philosophy, the theory is developed to identify the 
transition states and to classify them both dynamically and physically; see [13j [15] 
from application point of view of the theory. 

With this theory, many longstanding phase transition problems are either solved 
or become more accessible. The modeling and the analysis with the applications 
of the theory to specific physical problems, on the one hand, provide verifications 
of existing experimental and theoretical studies, and, on the other hand, lead to 
various new physical predictions. For example, our study [23 E] of phase 
transitions of both liquid helium-3 and helium-4 leads not only to a theoretical 
understanding of the phase transitions to superfluidity observed by experiments, 
but also to such physical predictions as the existence of a new superfluid phase C for 
liquid helium-3. Although these predictions need yet to be verified experimentally, 
they certainly offer new insights to both theoretical and experimental studies for a 
better understanding of the underlying physical problems. 
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2. Atmospheric Model over the Tropics 

Upwelling and zonal circulation over the tropics contains six convective cells, as 
shown in Figure |2.1| Usually the Walker circulation is referred to the Pacific cell 
over the equator, responsible for creating ocean upwelling off the coasts of Peru and 
Ecuador. It was discovered by Jacob Bjerknes in 1969 and was named after the 
English physicist Gilbert Walker, an early-20th-century director of British obser- 
vatories in Indian, who discovered an indisputable link between periodic pressure 
variations in the Indian Ocean and the Pacific, which he termed the "Southern 
Oscillation" . 
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Figure 2.1. Global atmospheric Walker circulation over the tropics. 



The Walker cell is of such importance that its behavior is a key factor giving 
rise to the El Nino (more precisely, the El Nino-Southern Oscillation (ENSO)) 
phenomenon. When the convective activity weakens or reverses, an El Nino phase 
takes place, causing the ocean surface to be warmer than average, reducing or 
terminating the upwelling of cold water. A particularly strong Walker circulation 
causes a La Nina event, resulting in cooler ocean temperature due to stronger 
upwelling; see Figure |2.2| 
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Figure 2.2. A schematic diagram of the Walker circulation at the 
equatorial Pacific Ocean, a main factor of the El Nino and the la 
Nina phenomena. 



The atmospheric motion associated with the Walker circulation affects the loops 
on either side. Under normal years, the weather behaves as expected. But every few 
years, the winters becomes unusually warm or unusually cold, or the frequency of 
hurricanes increases or decreases. This entirely ocean-based cell is seen at the lower 
surface as easterly trade winds which drives the seawater and air warmed by sun 
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moving towards the west. The western side of the equatorial Pacific is characterized 
by warm, wet low pressure weather, and the eastern side is characterized by cool, 
dry high pressure weather. The ocean is about 60 cm higher in the western Pacific 
as the result of this circulation. The air is returned at the upper surface to the 
east, and it now becomes much cooler and drier. An El Nino phase is characterized 
by a breakdown of this water and air cycle, resulting in relatively warm water and 
moist air in the eastern Pacific. Meanwhile, a La Nina phase is characterized by 
the strengthening of this cycle, resulting in much cooler water and drier air than 
the normal phase. During an El Nino or La Nina event, which occurs about every 
3-6 years, the weather sets in for an indeterminate period. 

It is then clear that the ENSO is a coupled atmospheric-ocean phenomena. The 
main objective of this section is to address the mechanism of atmospheric Walker 
circulations over the tropics from the dynamic transition point of view. 

The atmospheric motion equations over the tropics are the Boussinesq equations 
restricted on 9 = 0, where the meridional velocity component ug is set to zero, and 
the effect of the turbulent friction is taking into considering using the scaling law 
derived in [20] : 

du v _ u v u z ( . 2 du z 2u v 



(2.1) 



— — _ -(u ■ V)u v h v Au v -r 

at a \ a z dip 

on 1 d P 

p a dip 

du z „, u 2 ( . 2 du^ 2u 

- -(u- V)u z + -£ + v[Au. 



dt a \ a 2 dip a 2 

1 dp 
Po dz 



a 1 u z + 2Q,u v - — — - (1 - a T (T - T ))g, 



dT 

— = -(u-V)T+kAT, 

1 du v duz_ _ Q 
a dip dz 

Here Oi — Cih 2 (i = 0,1) represent the turbulent friction, a is the radius of the 
earth, the space domain is taken as M — x (a, a + h) with being the one- 
dimensional circle with radius a, and 

. „, u v d d . 1 d 2 d 2 

a dip dz a z d<p l oz 1 

For simplicity, we denote 

{x\,x?) = (cup, z), (ui,u 2 ) = (u v ,u z ). 

In atmospheric physics, the temperature T\ at the tropopause z = a + h is a 
constant. We take To as the average on the lower surface z = a. To make the 
nondimcnsional form, let 

x — hx' , t = h 2 t' / k, u = kv! /h, 

T = (T - T^T' + To - (T - T )x' 2 , 

p = Kvp p'/h 2 - gp Q (hx' 2 + a T (T - Ti)hx' 2 /2). 
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Also, we define the Raylcigh number, the Prandtl number and the scaling laws by 

, 00 % d argiTo-T^h 3 v . 

(2.2) R = , Pr = , dj = dli jv (i = 0, 1). 



Omitting the primes, the nondimensional form of (2.1 1 reads 
dm 



(2.3) 



dt 



du 2 
~dt 



Pr 



2 du 2 2 

Ami H 77 ui - b ui 

r Q dxi r 



dp 



LUU 2 ~ (U • V)«l 



1 



-UiU 2 , 



Pr 



2 dm 2 

Au 2 u 2 — biu 2 + RT 

r dxi r 

1 



dp 
dx 2 



+ Will — (u ■ V)«2 Ui , 

9xi 



= AT + u 2 - (u ■ V)T, 
= 0, 



where (xi, x 2 ) £ M = (0, 27rro) x (ro, ro + 1), <5o an d <5i are as in (2.2 1, (it • V) and 
A as usual differential operators, and 

2ft/i 2 

(2.4) uj= . 

K 

The problem is supplemented with the natural periodic boundary condition in 
the a; i -direction, and the free-slip boundary condition on the top and bottom bound- 
ary: 

(2.5) (u,T)(x 1 +2Trr ,x 2 ,t) = (u,T)(x 1 ,x 2 ,t), 



(2.6) 



«2 = 0, 77 — = 0, T = ip(xi) 
ax 2 

u 2 - 0, ^ = 0, T = 
dx 2 



at x 2 = ro, 
at x 2 = r + 1 . 



Here <p(x\) is the temperature deviation from the average To on the equatorial 
surface and is periodic, i.e., 

<p(xi)dxi — and <p(xi) = ip(x\ + 27rro). 

The deviation tp(x\) is mainly caused by a difference in the specific heat capacities 
between the sea water and land. 



3. Walker Circulation under the Idealized Conditions 



In an idealized case, the temperature deviation ip vanishes. In this case, the 
study of transition of (2.3 1 is of special importance to understand the longitudinal 
circulation. Here, we are devoted to discuss the dynamic bifurcation of (2.3), the 
Walker cell structure of bifurcated solutions, and the convection scale under the 
idealized boundary condition 



(3.1) 



ip(xi) = for any < x% < 2irr . 
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The following theorem provides a theoretical basis to understand the equatorial 
Walker circulation. The proof of this theorem follows the same method as given in 
where similar results are obtained for the classical Bernard convection. 



The orem 3.1. Under the idealized condition (3.1), the problem (2.3) with (2.5) 
and (2.6) undergoes a Type-I transition^ at the critical Rayleigh number R = R c . 
More precisely, the following statements hold true: 

(1) When the Rayleigh number R < R c , the equilibrium solution (u, T) = is 
globally stable in the phase space H defined by: 

H={(u,T) G L 2 (M) 3 | divu = 0, u 2 \ X2=ro ,r +i = 0, (u, T) satisfies (g}, 

where L 2 (M) stands for the square-integrable functions in M. 

(2) When R c < R < R c + 5 for some S > 0, this problem bifurcates from 
((u,T), R) = (0,i? c ) to an attractor Ar = S , consisting of steady state 
solutions, which attracts H\T , where T is the stable manifold of (u,T) = 
with codimension two. 

(3) For each steady state solution ipn = (ur,Tr) £ TZ,ur is topologically equiv- 
alent to the structure as shown in Figure \3~1\ 

(4) For any initial value tpo — (wo,Io) £ H \ (T U E), there exists a time 
to > such that for any t > to the velocity field u(t,ipo) is topologically 



equivalent to the structure as shown in Figure 3.2 either (a) or (b), where 
Tp = (u(t,4 } o),T(t,tpo)) is the solutions of the problem with V'(O) = ipo, and 

rro+l 

E = {{u,T)eH\ / Ul dx 2 = 0}. 




Figure 3.1. The cell structure of the steady state solutions in the 
bifurcated attractor Ar. 




(a) (b) 



Figure 3.2. The Walker cell structure with the cells separated by 
a cross channel flow: (a) a westbound cross channel flow, and (b) 
an eastbound cross channel flow. This cross-channel flow pattern 
has the same topological structure as the Walker circulation over 
the tropics and the Branstator-Kushnir waves in the atmospheric 
dynamics [H[5]. 



See |13l for precise definition. 
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A few remarks are now in order: 

First, in a more realized case where T = Tg + ip(xi) at xi = fo, the temperature 
deviation <p(x\) is small in comparison with the average temperature gradient To — 
T\. Therefore, the more realistic case can be considered as a perturbation of the 
idealized case. 

Second, mathematically, the conclusion that Ar — S 1 consisting of steady state 
solutions is derived from the invariance of (2.3 1 under a translation X\ — > X\ + a. 
It implies that for any two solutions ipi,ip2 £ Ar, they are the same up to a phase 
angle a in the longitude, namely 

ip2(xi,X2) — i/ji(xi +a,X2), for some a £ [0, 27iTo]. 

For the idealized case, this conclusion is nature as the equator is treated as homo- 
geneous. 

In addition, Assertion (2) amounts to saying that as the initial value V'o varies, 
i.e., an external perturbation is applied, the rolls of the Walker circulation will be 
translated by a phase angle a. This behavior is termed as the translation oscillation, 
which can be used to explain the ENSO phenomenon for the idealized circumstance. 

Third, when a deviation f(x\) is present, the translation symmetry is broken, 
and consequently, the bifurcated attractor Ar consists of isolated equilibrium solu- 
tions instead of a circle S 1 . Thus, the mechanism of the ENSO will be explained as 
state exchanges between these equilibrium solutions in Ar, which are metastable. 
We shall address this point in detail later. 

Fourth, by the structural stability theorems in [TO] , we see from Assertion (4) 
that the roll structures illustrated by Figure 3.2 (a) and (b) are structurally sta- 
ble. Hence, for a deviation perturbation ip(xi), these structures are not destroyed, 
providing a realistic characterization of the Walker circulation. 

Fifth, in Assertion (3), these rolls in the bifurcated solutions are closed. How- 
ever, under some external perturbation, a cross channel flow appears which sepa- 
rates the rolls apart, and globally transports heat between them. 

Sixth, it is then classical [TTJ [3J H] to show that the critical Rayleigh number 
is given by 

(a 2 + n 2 ) 3 + TT 2 {a 2 + TT 2 )S ~ 



(3.2) R c 



{a 1 + tH)(5i + 



Let g(x) = (x + tt 2 )5i + {{x + tt 2 ) 3 + ir 2 (x + ir 2 )5 )/x. Then, 
g'(x) — 0. Thus we get 

(3.3) S ia 4 c ~ (a 2 c + tt 2 ) 3 - Tt 2 (a 2 c + tt 2 )S + 'Sa 2 c (a 2 c + tt 2 ) 2 + it 2 a 2 c 5 a = 0. 

We assume that £i > Sq for h ^ 0. By a 2 — ir 2 /L 2 , when h is large, <5i 3> 1 and 
5\ ^> Sq, from (3.2) and (3.3) we derive that 

(3.4) a c = , 

Ol 

(3.5) L 2 c ^sY 2 /(n 2 + S Q y/ 2 , 

(3.6) Rc^^Si+^SoLl, 



These formulas ( 3.4 )-( 3.6 ) are very useful in studying large-scale convection motion, 
as shown in the next section. 
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4. Walker circulation under natural conditions 



4.1. Physical parameters and effects of the turbulent friction terms. It 

is known that air properties vary with the temperature and pressure. The table 
below lists some common properties of air: density, kinematic viscosity, thermal 
diffusivity, expansion coefficient and the Prandtl number for temperatures between- 
100°C and 100°C. 



Temperature 
t 


Density 
P 

(kg/m 3 ) 


Kinematic 
viscosity 

V 

(m 2 /s) X 10~ 6 


Thermal 
Diffusivity 

K 

(m 2 /s) X 10" 6 


Expansion 
coefficient 
a 

(1/k) X 10~ 3 


Prandtl 
number 
Pr 


-100 


1.980 


5.95 


8.4 


5.82 


0.74 


-50 


1.534 


9.55 


13.17 


4.51 


0.725 





1.293 


13.30 


18.60 


3.67 


0.715 


20 


1.205 


15.11 


21.19 


3.40 


0.713 


40 


1.127 


16.97 


23.87 


3.20 


0.711 


60 


1.067 


18.90 


26.66 


3.00 


0.709 


80 


1.000 


20.94 


29.58 


2.83 


0.708 


100 


0.946 


23.06 


32.8 


2.68 


0.703 



The average temperature over the equator is about 20 ~ 30° C. Based on the data 
in Table |4.1[ we take a set of physical parameters of air for the Walker circulation 
as follows: 

v = 1.6 x lCT 5 m 2 /s, k= 2.25 x lCT 5 m 2 /s, 
a T = 3.3 x lCT 3 / ^ p r = 0.71. 
The height h of the troposphere is taken by 

h = 8x 10 3 m. 



Thus, the Raylcigh number for the equatorial atmosphere is 
(4.1) R 



gaM-T,) ^ = 3 g x 1Ql9 _ a 



{a 2 + it 2 

(X 2 



= = 675, 



In the classical theory, we have the critical Rayleigh number 

(4.2) i? r =mm 
and the convective scale 

(4.3) L 2 = n 2 /a 2 c = 2. 
Namely, the diameter of convective roll is 

d = L c x h = 20km. 

However, based on atmospheric observations, there are six cells overall tropics, and 
it shows that the real convective scale is about 



(4.4) 



L = — = 6600 km . 



As a comparison, the values (4.2 1 and (4.3) from the classical theory are too 



small to match the realistic data given in (4.1 1 and (4.4 1 



We use ( 3.4 )-( 3.6 1 to discuss this problem. With the values in (4.1) and (4.4 1 
we take 



(4.5) 



S 1 = 2.7 x 10 2 



5 Q = 3.5 x 10 s 
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In comparison with (2.2), for air, the constants Co and C\ take the following values 

C = 1.37 x 1CT 12 rrT 2 ■ s" 1 , C x = 1.05 m~ 2 ■ s" 1 . 
Then, the convective scale L c in ( |3.5| takes the value 

(4.6) T c = [<5i/(5 ] 1/4 /i = 7500km, 

which is slightly bigger than the realistic value L = 6300km. And, the critical 
Rayleigh number in (3.6 1 takes 

R c tt 2 S 1 + tt 2 S L 2 c = 2.7 x 10 21 . 

Thus, the critical temperature difference is 

(4.7) AT C = 60°C, 

which is realistic, although it is slightly smaller than the realistic temperature 
difference on the equator. 

However, with the classical theory with no friction terms, the critical values 
would be L c = 20km and AT C = f 0~ 17o C, which are certainly unrealistic. 

4.2. Transition in natural conditions. We now return the natural boundary 
condition 

In this case, equations (2.3) admits a steady state solution 

(4.8) i> = (V,J), (V=(Vi,V a )). 
Consider the deviation from this basic state: 

u^u + V, T -^T+J. 



Then (2.3 1 becomes 

dUl , t Y7\ , UlU2 V [A x i d P I 

— - + (u ■ V)in H = Pr [Am - d ui - - — 

at r ox\ 

Vi V 2 
~ (u ■ V)ui — (u ■ V)Vi U2 Ui, 

(4.9) ^(..V^-^Pr^-WKT-jL] 

- (u ■ V)u 2 - (u ■ V)V 2 H L ui, 

ro 

dT 

— + (u ■ V)T = AT + u 2 - (u ■ V) J - (u ■ V)T, 
divu = 0. 

The boundary conditions are the free-free boundary conditions given by 
(u,T) is periodic in x\ — direction, 

(4.10) dm 

T = 0, u 2 = 0, - — = at x 2 = r ,r + I. 
ax 2 

It is known that |</?(xi)|, with AT = To — T\ = 100°C as unit, is small. Hence, 
the steady state solution (V, J) is also small: 

\\(V,J)\\v, =£< 1. 



Thus, (4.9 1 is a perturbated equation of (12. 3b. 
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Since perturbation terms involving (V, J) are not invariant under the zonal trans- 
lation (in the ^-direction), for general small functions <p(xi) 7^ 0, the first eigen- 



values of (4.9) are (real or complex) simple, and by the perturbation theorems in 



[H], all eigenvalues of linearized equation of (4.9 1 satisfy the following principle of 
exchange of stability (PES): 

for R < R s c , 

Re(3 e (R) { = for R = R e c , for any 1 < i < m, 

for R> R e c , 

for any j > m + 1 , 

where m = 1 as (3f(R) is real, m = 2 as (3f(R) is complex near R e c , and R e c is the 




critical Rayleigh number of perturbed system (4.9 1. 

The following two theorems follow directly from Theorem |3 . 1 1 and the perturba- 
tion theorems in [5]. 

Theorem 4.1. Let 0f(R) near R = R 6 C be a real eigenvalue. Then the system 



(4-9) has a transition at R = R £ c , which is either mixed (Type-Ill) or continuous 
(Type-I), depending on the temperature deviation ip{x{). Moreover, we have the 
following assertions: 

(1) // the transition is Type-I, then as R s c < R < R e c + 6 for some 6 > 0, the 
system bifurcates at R £ c to exactly two steady state solutions tpi and ip2 in 
H , which are attractors. In particular, space H can be decomposed into two 
open sets U\,U2'- 

H = U 1 + U 2l U l nU 2 = <j), -4> = 0edUiOdU 2 , 

such that ipi G Ui (i = 1,2), and ipi attracts Ui. 

(2) // the transition is Type-Ill, then there is a saddle-node bifurcation at R = 
R* with R* < R s c such that the following statements hold true: 

(a) if R* < R < R e c + 8 with R ^ R e c , the system has two steady state 
solutions ip 1 ? and ip^: which are attractors, as shown in Figure 
such that 

V>+ =0 forR*<R< R £ c . 

(b) There is an open set U C H with € U which can be decomposed into 
two disjoint open sets U = US; + U^} with tp± <= U± and ip± attracts 

Ug. 

(3) For any initial value ipo = (uq, To) £ for R > R* or ipo = (uq, Tq) S U? 
for R > R e c , there exists a time t^>Q such that for any t > to the velocity 
field u(t, ipo) is topologically equivalent to the structure as shown in Figure 



5.3 



3.2 either (a) or (b), where ip = (u(t,tpo),T(t,rpo)) is the solutions of the 



problem with ip(0) = tpo- 

Theorem 4.2. Let Pl(R) be complex near R = R e c . Then the system {^.9\ bifur- 
cates from (ip,R) = (0, R 5 .) to a periodic solution ipii(t) on R e c < R, which is an 
attractor, and ipn(t) can be expressed as 

(4.11) ip R (t) = Anicospttpt + sinpi^i) + o(\A R \,e), 

where Ar = a(R— R e c ),oi and p are constants depending on tp(xx), and ip\,ipi are 



first eigen] 'unctions of linearized equations of (2.3) 
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A few remarks are now in order. 



First, Theorems |3.1| [4~T] and |4~2| provide the possible dynamical behaviors for 
the zonal atmospheric circulation over the tropics. Theorem |3.1| describes the trans- 
lation oscillation, and does not represent a realistic explanation to the ENSO. 



Second, the periodic solution (4.11 1 characterizes a roll pattern translating with 



a constant velocity over the equator, eastward or westward, as shown in Figure ET 



(a) and (b) The time-periodic oscillation obtained in Theorem 4.2 does not represent 
the typical oscillation in the ENSO phenomena, as the Walker circulation does not 
obey this zonal translational oscillation. 





(a) 



(b) 



Figure 4.1. Time-periodic translation of the Walker circulation 
pattern; (a) an eastward translation, (b) a westward translation. 



Third, Theorem |4 . 1 1 characterizes the oscillation between metastable states, The 
oscillation between the metastable states in Theorem |4.1| can be either Type-I or 
Type-Ill, depending on the number 



where is the first eigenvector of the linearized equation of (4.9 1 at R = R 



*Sf* is the corresponding first eigenvector of the adjoint linearized problem, and G 



represents the nonlinear terms in (4.9 1 



Namely, if b 0, the transition is Type-Ill and if b = 0, the transition is Type- 
I. From mathematical viewpoint, for almost all functions ip, b =/= 0. Hence, the 
case where 6^0, consequently the Type-Ill transition in Theorem |4.1[ is generic. 
In other words, the Type-Ill transition derived in this theorem provides a correct 
oscillation mechanism of the ENSO between two metastable El Nino and La Nina 
events, and we shall explore this point of view in the next section in detail. 

5. Metastable Oscillation Theory of ENSO 

It is well known that the Walker cell at the equatorial Pacific Ocean is closely 
related to the ENSO phenomenon. Its behavior is the key to the understanding of 
ENSO. 
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Figure 5.1. SOI indicates 
the state of the Walker 
circulation. 



Southern Oscillation Phenomenon. The Southern Oscillation Index (SOI) 
gives a simple measure of the strength and phase of the Southern Oscillation, and 
indicates the state of the Walker circulation. The SOI is defined by the pressure 
difference between Tahiti and Darwin. When the Walker circulation enters its El 
Nino phase, the SOI is strongly negative; when it is in the La Nina phase, the SOI 
is strongly positive; and in the normal state the SOI is small; see Figure |5.1| 

This point can be further examined by further observational data of SOI as 



shown in Figure 5.2 In the SOI diagram, we observe that there are two groups 
of oscillations: the relatively large amplitude fluctuation and the relatively small 
amplitude fluctuation. The large one occurred in 1950-1951, 1955-1956, and 1974- 
1975 for positive SOI higher than 13, in 1965-1966, 1972-1973, 1977-1978, 1981- 
1982, and 1987-1988 for negative SOI lower than -12. The small one with SOI 
between-8 to 8 appeared in other years. 



SOUTHERN OSCILLATION INDEX 

1950 to 1999 
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Year 19xx (Monthly Data) 



Figure 5.2. SOI diagram, adopted from Web: BOM-Monitoring 
Climate-Climate variability and El Nino. 



El Nino and La Nina States. The non-homogeneous temperature distribu- 
tion shows that the locations of the cells over the tropics are relatively fixed. This 
suggests that the homogeneous case as described by Theorem |3.1| is not a valid de- 
scription of the realistic situation. In addition, the translation oscillation formation, 
suggest in Theorem |4.2| does not describe the ENSO either. 

Hence the nature theory for the atmospheric circulation over the tropics asso- 



ciated with the Walker circulation and the ENSO is given by Theorem 4.1 It is 



known that the steady state solutions of (2.3 1 can be written as 

(5.1) 4> + i> R , 
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where %j) — (V, J) is the steady state solution given by (4.8), and ip^, are the stable 
equilibria of (14.91), derived in Theorem 



In fact, the ENSO phenomenon can 



4.1 



be explained as the transition between the two metastable steady state solutions 



Theorem 4.1 amounts to saying that for a mixed transition, there are two critical 
Rayleigh numbers R* and R% with R* < R% where R* = R c = 2.7 x 10 21 . At 
R = R* , the equation (4.9 1 has a saddle-node bifurcation, and at R — R E C has the 



normal transition, as shown in Fig ure |5.3| When the Rayleigh number R is less 
than R* , i.e. R < R *, the system ( |2.3[ ) has only one stable steady state solution 
tp = (V, J) as in (4.8 1 which attracts H. When R* < R < R e c , this system has two 
stable equilibriums 

(5.2) 



and tp + 6 Uj?, for R* < R < R £ c 



where U+, U 1 } C H are basins of attraction of ip and ip + ip R . And when R e c < R, 
this system has also two stable equilibriums 



(5.3) 



ip + ipne U R and $ + ip± € U? for R* c < R, 



with U R and are their basins of attraction. 




El Nino phase 



R Rayleign number 



i[>" La Nina phase 



Figure 5.3. The transition diagram for the equatorial motion 



equations (2.3 1, R* is the saddle-node bifurcation point, R e c is the 



first critical Rayleigh number. 



Since the problem (2.3 1 with a natural boundary condition is a perturbation for 



the idealized condition, ip — {V, J) is small, and the velocity field V is almost zero. 
Hence the equilibrium 



ip for R<R* C , 

_^ + ^+ for R>R* C , 
represents the El Nino phase, and the equilibrium 



represents the normal phase for R* < R < R* c , and the La Nina phase for R* < R. 
It is clear tl 
respectively. 



It is clear that both are metastable, with as their basins of attractions 



Oscillation Mechanism of ENSO. The above theoretical studies suggest that 
ENSO is the interplay between two oscillation processes. The first is the oscillation 
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between the metastable warm event (El Nino phase, represented by and cold 
event (La Nina phase, represented by ^f^). The second oscillation is the oscillation 
of the Rayleigh number caused essentially by the some spatiotemporal oscillation 
of the sea surface temperature (SST) field. 

Here we present a brief schematic description on the interplay between these two 
oscillation processes of the ENSO, based on the saddle-node transition diagram in 



Figure 5.3 rigorously proved in Theorem 4.1 



We start with three physical conclusions: 

1 We observe that as R decreases, the normal and the La Nina phase 
weakens and its basin of attraction shrinks (to zero as R approaches to 
R*). 

2 As R increases in the interval R% < R, the strength of the El Nino phase 
increases. As R increases in (i?*,i?J), the basin of attraction of the El 
Nino phase shrinks. In particular, near R e c , the El Nino phase is close to 
the interaction of the two basins of attraction of the El Nino and La Nina 
phase; consequently, forcing the transition from the El Nino phase to the 
La Nina phase. 

3 Also, we see that as the El Nino event intensifies, the SST increases, 
and as the normal and La Nina event intensifies, the SST decreases. 

From the above three physical conclusions, we obtain the following new mecha- 
nism of the ENSO oscillation process: 

I When R* < R < R e cl represents the normal condition, and the corre- 
sponding upwelling near Peru leads to the decreasing of the SST, and leads 
to R approaching R*. On the other hand, near R*, the basin of attraction 
of the normal condition shrinks, and due to the uncertainty of the ini- 
tial data, the system can undergo a dynamic transition near R* toward to 
the El Nino phase 

II For the El Nino phase = ip near R* , however, the velocity V is small, 
and SST will increases due to the solar heating, which can not be trans- 
ported away by weak oceanic currents without upwelling. Hence the corre- 
sponding Rayleigh number R increases, and the El Nino phase intensifies. 
When the Rayleigh number approaches to the critical value R e cl with high 
probability, the system undergoes a metastable transition from the El Nino 
phase to either the normal phase or the La Nina phase depending on the 
Rayleigh number R and the the strength of the phase 
Ill With delaying effect, strong El Nino and La Nina occur at the Rayleigh 
number R larger than the critical number R e c , and the system will repeat 
the process described in items I-II above. 

In summary, this new mechanism of ENSO as an interplay of the two processes 
leads to both the random and deterministic features of the ENSO, and defines a new 
natural feedback mechanism, which drives the sporadic oscillation of the ENSO. 
The randomness is closely related to the uncertainty /fluctuations of the initial data 
between the narrow basins of attractions of the corresponding metastable events, 
and the deterministic feature is represented by a deterministic coupled atmospheric 
and oceanic model predicting the basins of attraction and the SST. It is hoped 
this mechanism based on a rigorous mathematical theory could lead to a better 
understanding and prediction of the ENSO phenomena. 
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